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yielding positive-definite and monotonically decreasing gluon masses, in agreement with a variety 
of less formal considerations. 
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1. Introduction 

It is well-established by now that the dynamical generation of an effective gluon mass ^ ex- 
plains in a natural and self-consistent way the infrared finiteness of the (Landau gauge) gluon prop- 
agator and ghost dressing function, observed in large- volume lattice simulations for both SU (2) 
and SU (3) gauge groups Given the nonperturbative nature of the mass generation mechanism, 
the Schwinger-Dyson equations (SDEs) constitute the most natural framework for studying such 
phenomenon in the continuum [Q, ^, ^, Specifically, we will work in the framework provided by 
the synthesis of the pinch technique (PT) [|l|, [Sj, with the background field method (BFM) [0], 
known in the literature as the PT-BFM scheme [^, [IT]]. 

Probably the most crucial theoretical ingredient for obtaining out of the SDEs an infrared-finite 
gluon propagator, without interfering with the gauge invariance of the theory, encoded in the BRST 
symmetry, is the existence of a set of special vertices, to be generically denoted by V and called 
pole vertices. These vertices contain massless, longitudinally coupled poles, and must be added 
to the usual (fully dressed) vertices of the theory. They capture the underlying mass generation 
mechanism, which is none other than a non-Abelian realization of the Schwinger mechanism. In 
addition to triggering the Schwinger mechanism, the massless poles contained in the pole vertices 
act as composite, longitudinally coupled Nambu-Goldstone bosons, maintaining gauge invariance 
and preserving the form of the Ward identities (Wis) and the Slavnov-Taylor identities (STIs) of the 
theory in the presence of a dynamically generated gluon mass. In fact, recent studies indicate that 
the QCD dynamics can indeed generate longitudinally coupled composite (bound-state) massless 



poles, which subsequently give raise to the required vertices V []12|, |13|]. 

At the level of the SDEs, the analysis finally boils down to the derivation of an integral equa- 
tion, to be referred as the mass equation, that governs the evolution of the dynamical gluon mass, 
m^(q^), as a function of the momentum q^. The main purpose of this presentation is to report on 
recent work [jl^, where the complete mass equation has been obtained in the Landau gauge em- 
ploying the full SDE of the gluon propagator and using as a guiding principle the special properties 



of the aforementioned vertices V (for related studies in the Coulomb gauge see, e.g., [|15|, |160). 
In this context, the detailed numerical solution of the fuU mass equation (for arbitrary values of 
the physical momentum), reveals the existence of positive-definite and monotonically decreasing 
solutions. 

2. The SDE of the gluon propagator 

The full gluon propagator A''^y{q) = 5"''A^y{q) in the Landau gauge is given by the expression 

A^y{q) = -iP^y{q)A{q^y, P^^(q) = g^^ (2.1) 

and its inverse gluon dressing function, J{q^), is defined as 

A-\q^)=q'j{q'). (2.2) 

The usual starting point of our dynamical analysis is the SDE governing the gluon propagator. 
Specifically, within the PT-BFM formalism, one can consider the propagator connecting a quantum 
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Figure 1: The SDE obeyed by the QB gluon propagator. Black blobs represents fully dressed 1-PI vertices; 
the small gray circles appearing on the external legs (entering from the right, only!) are used to indicate 
background gluons. 



(Q) with a background (B) gluon, to be referred as the QB propagator and denoted by A{q-^). The 
SDE of the above propagator is shown in Fig. |l|, and it may be related to the conventional QQ 
propagator, A(^^), connecting two quantum gluons, through the powerful background-quantum 
identity [|, [nil 

A{q^) = [l+G{q')]A{q^), (2.3) 

In this identity, the function G{q^) corresponds to the g^v form factor of a well known two-point 
function [Q, ^ [T^. Then, the corresponding version of the SDE for the conventional gluon propa- 
gator (in the Landau gauge) reads [| , 11 ] 



A-'(r)P,A^ = '^''"'^f^^f'^" . (2.4) 

where the diagrams (a;) are shown in Fig. [!} The relevant point to recognize here is that the 
transversality of the gluon self-energy is realized according to the pattern highlighted by the boxes 
of Fig. [T], namely, 

/[(ai) + (a2)]^v=0; ^^[(as) + (a4)]^v = 0; ^^[(aj) + K)]mv = 0. (2.5) 



3. Derivation of the gluon mass equation 

As has been explained in detail in the recent literature [|^, 13], the Schwinger mechanism 
allows for the emergence of massive solutions out of the SDE, preserving, at the same time, the 
gauge invariance intact. At this level, the triggering of this mechanism proceeds through the in- 



clusion of the pole vertices V in the SDE Eq. ( |2.4| ). From the kinematic point of view, we will 



describe the transition from a massless to a massive gluon propagator by carrying out the replace- 
ment (Minkowski space) 



A-\q')=q^J{q^) A;„\q^) = q^Uq^) -m\q'). (3.1) 
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Notice that the subscript "m" indicates that effectively one has now a mass inside the correspond- 
ing expressions: for example, whereas perturbatively J{q^) ~ In^^, after dynamical gluon mass 
generation has taken place, one has Jm{q^) ~ ln(^^ + m^). Then, gauge invariance requires that the 
replacement given in Eq. (3J_) be accompanied by the following simultaneous replacement of all 
relevant vertices 

T^Y' = Y,n + V, (3.2) 

where V must be such that the new vertex Y' satisfies the same formal Wis (or STIs) as Y before. 
The most familiar case is that of the BQ^ vertex F'^^y, whose pole part must satisfy the WI 

q''Vanv{q,r,p)=m^{r')P^y{r)-m'{p')P^^{p), (3.3) 

when contracted with respect to the momentum of the background gluon. In complete analogy with 
the above case, one may use the WI satisfied by the conventional BQ^ vertex, namely, 

+ r'7^'%vp(r,p,<? + 0], (3.4) 
in order to deduce that, after the replacement Eq. (^^, its V part satisfies [ p^ ] 

q^'Ktvpiq.W) = ig'[f'''r'Vyp,{p,t,q + r)+nf'">Vp^,it,r,q + p) 

+ r'^r'%vp(r,p,<?+0]- (3.5) 

Finally, as a large variety of lattice simulations and analytic studies suggest, we will take for granted 
that the ghost propagator D{q^) remains massless in the Landau gauge. The main implication of 
this property for the case at hand is that the (fully-dressed) BFM gluon-ghost vertex, appearing in 
graph (as), does not need to be modified by the presence of V-type vertices. 

Quite remarkably, the above Wis, supplemented by the totally longitudinal nature of the pole 
vertices, are the only properties that one needs for deriving the mass equation; in particular, the 
closed form of the pole vertices is not needed. 

According to the previous discussion, after the inclusion of the pole vertices, the gluon SDE 
Eq. ( lA ) becomes in the Landau gauge 

[q J^{q)-m [q )\P^y{q) = l+G{q^) ' 



where the prime indicates that (in general) one must perform the simultaneous replacements Eq. (3.1) 
and Eq. (3.2) inside the corresponding diagrams. Evidently, the Ihs of Eq. (3.6) involves two un- 
known quantities, Jm{q^) and m^{q^), which will eventually satisfy two separate, but coupled, 
integral equations of the generic type 



Jm{q^) = 1+ j^J(f\{q^,n?Am), 

) = J^,je2iq\m\A,„), (3.7) 



m^(q^ 



such that J^,^ / 0, as q^ — )• 0. In order to derive the closed form of the mass equation Eq. (3.7), 
one must identify all mass-related contributions coming from the V vertices that are contained 
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Figure 2: Diagrammatic representation of the condensed operations leading to the all-order gluon mass 
equation, where we have introduced the shorthand notation in^{q^) = m^{q^)[\ + G{q^)\. All internal prop- 
agators are in the Landau gauge. 



in the Feynman graphs comprising the rhs of Eq. (3^). With the transversality of both sides of 



Eq. ( |3.6[ ) guaranteed by the presence of the pole vertices, it is far more economical to derive the 
mass equation by isolating the appropriate cof actors of q^qy/q^, to be denoted by of (q^), on both 



sides. Notice that selecting the g^y, or taking the trace in Eq. ( |3.6[ ), would entail the use of the 
special seagull identity [|^, 20]. 

The most important steps of this construction may be summarized as follows: (i) From the 
previous comments about the BFM gluon-ghost vertex, graph (cij) have not V-component. (ii) The 
WI Eq. (3.5 ) and the longitudinality condition for the BQ^ pole vertex may be used to demonstrate 
that the V-component of graph {a'^) vanishes in the Landau gauge, (iii) The contribution a^{q^) 
stems solely from the combination r^V of the product FT' appearing in graph (a'^). 

Thus, one concludes that the complete mass equation can be written as 



i[a\{q')+aW)] 
l+Giq^) 



(3.8) 



Interestingly enough, the entire procedure may be pictoriaUy summarized, in a rather concise way, 
as shown in Fig. ^ 



4. Complete mass equation and numerical results 



The final equation obtained from Eq. (3.8) reads (Euclidean space) 

.2. 



^^^^ ^-{qhsy-^qsqy) U\k^)C[Y{k + q)-Y{k)]Al{k)hy'{k + q), (4.1) 

Jk 



with C = ^hCaOCs, and 

Y{k^) = ^^^A«P(/)A/^-(/ + ^)r,p^, (4.2) 

corresponding to the subdiagram on the upper left comer of (ae) [see also Fig. ^]. 

Even though Eq. ( |4.1[ ) forms part of a system of coupled equations, see Eq. ( p.7| ), in what fol- 
lows we will study it in isolation, given that the corresponding equation for J,„{q^) is unknown. To 
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Figure 3: Left panel: Modification in the shape of the two-loop dressed kernel ^(y) when varying C. As C 
increases, the kernel effectively reverses its sign, showing a deep negative well in the low momenta region. 
Right panel: Typical monotonically decreasing solution of Eq. (5.1). The case shown has been obtained for 
the special values C — 0.88 and 1.85, corresponding value of the coupling as ~ 0.88 and 0.22, respectively. 
The solutions have been normalized so that at the origin they match the corresponding (Landau gauge) lattice 
value A-^O) ~ 0.141 GeV^, namely m(0) = 375 MeV. 



that end, we will treat the gluon propagators appearing in the mass equation as external quantities, 
using lattice results for their form [||]. 

In addition, the rhs of Eq. ( [4.2[ ) depends on the full three-gluon vertex r^yp^g , whose exact form 
is not known. We will therefore approximate Y{k^) by its one-loop expression, obtained by sub- 
stituting tree-level values for all quantities appearing in the integral; a lengthy but straightforward 
calculation yields [Euclidean space, momentum subtraction (MOM) scheme] 



i'sl'^') = -77^7 log ^- (4.3) 



15. k'^ 



(471)2 4 ^ 



After these considerations, and using spherical coordinates x = and y = k^, let us study the 



deep infrared limit a; — > of Eq. (4. 1 ), given by 



J(i{y) = {\\-lCY{y)\Z'{y)}'. ZM =yA(y). (44) 

Even though the value of C is fixed (see above), in what follows we will treat it as a free pa- 
rameter, in order to study what happens to the gluon mass equation when one varies independently 



as and C. The reason for doing this is that, whereas Eq. ( [^.3[ ) furnishes a concrete form for the two- 
loop dressed correction, by no means does it exhaust it; thus, by varying C, one basically tries to 
mimic further correction that may be added to the skeleton provided by the YR{k^) of Eq. ( p3| ) (for 
a fixed value of as); indeed rescaling C is equivalent to rescaling YR{k^). Of course, the real value 
of C will emerge as a special case of this general two-parameter study. Finally, it is convenient to 
define the "reduced" C, = C /JiUCa, and drop the suffix "r". 
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Let us first set C = 0, thus turning off tlie two-loop dressed contributions. Then, integrating 
( p!4l ) by parts, one obtains 

'"'(O) = o .""'wnM r dy[mHy)]'z\y) ■ (4.5) 



87r[l + G(0)]yo 

Given that 1 + G(0) > 0, it is clear that a monotonically decreasing gluon mass, namely [?m^(3')]' < 0, 
expected on physical grounds, would give rise to a negative value for m^(0), which is physically 
wrong. Thus, the only way to reconcile a positive-definite and monotonically decreasing gluon 
mass is to obtain an effective reversal of sign from the two-loop dressed contributions; as we will 
see, this is indeed what happens. 

To study this crucial point in detail, let us consider how the shape of the kernel changes 
as C is varied. In Fig. ^, one observes that, as C increases, ^ displays a less pronounced positive 
(respectively negative) peak in the small (respectively large) momenta region. Next, for C > 0.37, 
a small negative region starts to appear in the infrared, which rapidly becomes a deep negative 
well for y < 0.6, with becoming positive for higher momenta. Therefore we observe that the 
addition of the two-loop dressed contributions counteracts the effect of the overall minus sign of 
Eq. ( |4.4[ ), by effectively achieving a sign reversal of the kernel. Indeed, one concludes that there 
exists a critical value C « 0.56 such that, if C > C, Eq. ( p| ) will display at least one physical 
monotonically decreasing solution for a suitable value of the strong coupling as- 

Finally, to see if the picture sketched above is confirmed when x y^O, one can study numeri- 
cally the solutions of Eq. (4^1) following the algorithm described in [[l4l]. In this case the absence of 



solutions persists until the critical value C is reached, after which one finds exactly one monoton- 
ically decreasing solution. Specifically, in Fig. ^ we plot the solutions for the most representative 
C values. The value C = « 0.88 corresponds to the case in which Y is kept at its lowest order 
perturbative value, whereas C = 1.85 corresponds to the standard MOM value as = 0.22 at = 4.3 



GeV pip. As can be readily appreciated, the masses obtained display the basic qualitative features 
expected on general field-theoretic considerations and employed in numerous phenomenological 
studies; in particular, they are monotonically decreasing functions of the momentum and vanish 
rather rapidly in the ultraviolet. 

5. Conclusions 



In this presentation we have reported recent progress [14] on the study of the nonpertubative 
equation that governs the momentum evolution of the dynamically generated gluon mass. By ap- 
pealing to the existence of the special nonperturbative vertices V associated with the Schwinger 
mechanism, we have outlined the methodology that allows for a systematic and expeditious identi- 
fication of the parts of the SDE that contributes to the mass equation. The numerical analysis of the 
resulting mass equation reveals that the inclusion of two-loop dressed contributions has a profound 
impact on the nature of the mass equation, already at the qualitative level. Indeed, they are crucial 
in order to obtain physically meaningful solutions out of the mass equation, i.e., positive-definite 
and monotonically decreasing solutions for the effective gluon mass. 

In the future, given the importance of the term Y{k^) for this entire construction, it would be 
particularly important to determine its structure beyond the perturbative one-loop approximation 
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used. Nevertheless, even with the approximate version of Y{k^), the full mass equation provides a 
natural starting point for calculating reliably the effect that the inclusion of light quark flavors might 



have on the form of the gluon propagator, and complement recent studies based on the SDEs [ |22| ] 
as well as lattice simulations [|^]. 
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